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Abstract 

We study the statistical properties of eigenvalues of the Hessian matrix H 
(matrix of second derivatives of the potential energy) for a classical atomic 
liquid, and compare these properties with predictions for random matrix mod- 
els (RMM). The eigenvalue spectra (the Instantaneous Normal Mode or INM 
spectra) are evaluated numerically for configurations generated by molecular 
dynamics simulations. We find that distribution of spacings between nearest 
neighbor eigenvalues, s, obeys quite well the Wigner prediction s exp(—s 2 ), 
with the agreement being better for higher densities at fixed temperature. 
The deviations display a correlation with the number of localized eigenstates 
(normal modes) in the liquid; there are fewer localized states at higher den- 
sities which we quantify by calculating the participation ratios of the normal 
modes. We confirm this observation by calculating the spacing distribution 
for parts of the INM spectra with high participation ratios, obtaining greater 
conformity with the Wigner form. We also calculate the spectral rigidity and 

find a substantial dependence on the density of the liquid. 

PACS numbers: 75.10.Nr, 61.20.-p, 61.43.-j, 64.70.Pf, 05.40.-a, 05.45+b 
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I. INTRODUCTION: 



The local topography of the potential energy surface, as characterized by the ensemble 
averaged spectrum (INM spectrum) of eigenvalues of the second-derivative matrix (Hessian) 
of the potential energy function, have been studied in recent years as an approach to the 
analysis of dynamics in liquids . Broadly, the study of the INM spectra has been directed 
at analysing short time dynamics as in studying solvation 0], and at elucidating information 
about pathways to long time relaxation in the form of potential energy barriers etc [0,0 • 

Considerable effort has also been dedicated to developing analytical theories for calcu- 
lating the INM spectra within an equilibrium description |§-f|. The approach in much of 
these attempts has been to formulate the problem of calculating the INM spectrum as an 
exercise in random matrix theory. If one treats the individual elements of the Hessian ma- 
trix as independent, and distributed according to the appropriate Boltzmann weight, then 
the Hessian may be viewed as a real, symmetric random matrix with a known distribution 
of matrix elements. Two properties, however, distinguish the Hessian from the standard 
corresponding case treated in random matrix theory: (i) The diagonal entries of the Hes- 
sian are related to the off-diagonal entries by the property Tt°f = — J2j^i^tj > where i,j 
label the particles and a, [3 the spatial coordinates x, y, z. (ii) For liquids with short-ranged 
interaction potentials, the Hessian matrix is sparse, with the fraction of non-zero entries, p, 
depending on the system size as p ~ 1/N. 

In view of the above considerations, it is of interest to inquire to what extent the INM 
spectrum displays universal features idenitied in random matrix theory. In this paper, we 
address this question by obtaining INM spectra numerically for a model atomic liqiud that 
has been studied in the context of slow dynamics in supercooled liquids |f[0f| . The statistics 
we consider are the spacing statistics between nearest neighbor eigenvalues and the spectral 
rigidity, which we explain below. 
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II. EIGENVALUES OF THE HESSIAN 



The model liquid we study is a binary mixture || composed of 80% of particles of 
type A and 20% of type B, interacting via the Lennard- Jones potential, with Lennard- Jones 
parameters €ab/caa = 1-5, cbb/^aa = 0.5, ctab/caa = 0.8, and (JbbI&aa = 0.88, and a ratio 
of masses m^/m^ = 1. Lennard- Jones reduced units are used to report all the quantities, 
in terms of the A particle paramters €aa, &aa and m^: temperatures as T* = ksT/eAA, 

2 

densities as p* = pjo\ A and Hessian eigenvalues A* = \ mA<TAA . Further details may be 



found in Molecular dynamics simulations of the liquid are performed at ten reduced 
densities p* from p* = 0.65 to 1.40, at reduced temperature T* = 1.0. Two hundred sample 
configurations in each case are chosen from the equilibrated trajectory for the INM analysis. 
For each of these configurations, the Hessian is calculated and diagonalized numerically 
to obtain the eigenvalues \ as well as the eigenvectors ej. The eigenvectors are used to 
calculate the localization properties of the normal modes, via the participation ratio 
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The participation ratio thus defined is small (order of 1/N) for localized modes and large 
(order of 1) for extended modes. 

By constructing the histogram of eigenvalues A of the Hessian for all configurations 
considered, we obtain the INM Density of States (DOS) or spectrum [11]. Figure 1 shows 
the INM spectrum D(X) verses A for different densities. Note that the DOS is very different 
from the Wigner semi-circle distribution, obtained in RMM. We know from the literature 
|i~3| that the correlation and spacing functions are universal in certain regions no matter 
what the DOS is. Thus we use the corresponding statistics of spacings between eigenvalues, 
Wigner-Dyson statistics, as the standard of reference. 
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III. UNFOLDING THE SPECTRUM 



The statistical analysis of the numerical data proceeds by first using an unfolding proce- 
dure. The numerical calculation yields the eigenvalues of the Hessian which is ordered and 
forms the sample spectrum {Ai, A2, A„}. In order to analyse the spacing statistics, one 
must transform the eigenvalues Aj in such a way that the transformed eigenvalues Q are 
uniformly distributed. That is, the spectral density -D(C) — 1- This procedure is referred to 
as "unfolding" the spectrum fli~4l , |l"2"ll . In general such a transformation for a spectral density 
function -D(A) is most easily accomplished through its cumulative distribution 

C(A) = [ X D(\)d\, (III.l) 



by defining ((X) = C(A). For a discrete spectrum such as the ones we consider here, 
the corresponding procedure would be to extract from the sample spectrum the "smooth" 



part of the "staircase" cumulative distribution |]12| . In this work, the procedure we adopt 
to estimate the smooth part of the cumulative distribution is to evaluate the cumulative 
distribution for the union of all eigenvalues Aj obtained for 200 independent configurations 
at each density and temperature value. Such a smooth cumulative spectrum along with 
the "staircase" cumulative spectrum for a single configuration, is shown in Fig. 2, and the 
unfolding procedure is indicated. 



IV. SPACING DISTRIBUTION 

The spacing distribution P(s) for the random matrix models is defined as the probability 
of finding the next nearest neighbor eigenvalue of the spectrum to be at a distance s i.e. 
Si = A '+^~ A ' where A is the mean level spacing. In the present case where we use unfolded 
eigenvalues, Si = Ci+i — 0- Then P(s) = As^e~ Bs2 where (3 — 1 for the orthogonal random 
matrix model, which is the "Wigner surmise" [0. This spacing distribution arises in var- 
ious quantum systems which show an underlying classically chaotic behavior e.g. quantum 
billiards, quantum dots, nuclear spectra, disordered mesoscopic systems. The system we 
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analyze, on the other hand, is a classical liquid, with a disordered microscopic structure. 
The spacing distributions obtained are shown in Fig. 3 for three densities of the liquid. 
We find that to a first approximation, the cases considered display the universal behavior 
according to the Wigner surmise, with the agreement being better for the liquid at higher 
densities. We note here that in Ref. || a sparse random matrix is used to calculate the 
INM DOS analytically for a one- dimensional system which is very different from the usual 
Wigner semi-circle, and displays qualitative features very similar to the DOS we calculate 
numerically here. It would be very interesting to see whether the correlators of the sparse 
random matrix proposed in ref. || captures the above behaviour for different densities for 
the liquid. 



V. SPECTRAL RIGIDITY 

We next study the spectral rigidity, which measures the fluctuations of the number of 
eigenvalues in a window of given size as a function of the size of the window (or equivalently, 
the average number of eigenvalues expected in the window). The number fluctuations are 
plotted as a function of the mean number of eigenvalues in Fig. 4 for the same three densities 
as in the preceding figures. For the case of the Poisson spectrum and the harmonic oscillator 
14 the number fluctuations are < 5N 2 >~< iV > and < SN 2 >= constant, while for the 



Gaussian random matrix ensembles < 5N 2 >~ In < N > for N ^> 1. 

The number fluctuations are of the form < 5N 2 >~ iV 7 where 'j(p) depends on the 
liquid's density. This situation is reminescent of that found for the mobility edge in the 
quantum hall effect and in the Anderson metal-insulator transition. A detailed analysis of 
this awaits future work. 



VI. PARTICIPATION RATIO 



In order to get some insight into the reason for the systematically better agreement 
with increasing density of the spacing distribution, we consider the localization properties 
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of the normal modes in the liquid. In the standard random matrix case for orthogonal 
matrices, the eigenvectors are all extended, while from previous numerical and analytical 
studies we know that a fraction of the INM eigenmodes are localized. As described earlier, 
we use the participation ratio to quantify the localization of modes, averaging over modes 
corresponding to eigenvalues in each histogram bin, for the unfolded eigenvalues. In Fig. 5, 
the participation ratio is plotted as a function of the unfolded eigenvalues for T* = 1.00, for 
values of the density p* = 0.65, 1.00 and 1.35. One notes that for the highest density p = 1.35 
the participation ratios are highest overall, while for the lower densities the participation 
ratios are quite small for a substantial fraction of the eigenmodes, indicating a large number 
of localized modes. 

Next, we calculate the spacing distribution for p* = 1.35 for unfolded eigenvalues between 
0.1 and 0.6 for which the participation ratio is high (> 0.75) and relatively unchanging (Fig. 
|5[). The resulting spacing distribution is shown in Fig. 6, along with the spacing distribution 
for the entire eigenvalue spectrum and the expectation based on the RMM result. The 
data shown clearly demonstrate that the spacing distribution is practically identical to the 
standard RMM result, confirming the speculation that the increasing fraction of localized 
states at lower densities are responsible for the deviations at these densities from the standard 
RMM result. Indeed, this observation has been used, in a different context to locate the 



mobility edge in disordered systems [20 1 



VII. CONCLUSIONS 

We have presented the spacing statistics and spectral rigidity for numerically calculated 
INM spectra. The spacing statistics is seen to conform better with increasing density (at 
fixed temerature), with the predictions for random matrix models. We demonstrate that 
the source of deviations from RMM predictions is related to the presence of localized in- 
stantaneous normal modes in the liquid, whose number is greater for lower density. As the 
two features distinguishing the INM spectra from the standard random matrix case are the 
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nature of the diagonal elements and the sparseness of the Hessian matrix, further under- 
standing of the non-universality of the INM spectra are to be sought in the manner in which 
these aspects affect the INM spectral statistics. 
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FIGURES 

FIG. 1. INM Density of States shown for three densities. Inset shows the DOS on a logarithmic 
scale. Eigenvalues A are expressed in units of eA \ , and D(X) in units of A ° AA . 

FIG. 2. A portion of the staircase cumulative spectrum from a single configuration is shown, 
along with the estimate of the smooth part of the cumulative distribution, obtained from combining 
eigenvalues from 200 configurations, for density p* = 1.0 and temperature T* = 1.0. The arrows 
indicate the mapping of any given eigenvalue \ to the unfolded eigenvalue Q. The inset shows 
the cumulative spectra for the full range. Eigenvalues A are expressed in units of m e *Ji , and 
cumulative probabilities £ are dimensionless. 

FIG. 3. The level spacing distribution P(s) for densities p = 0.65, 1.0, 1.35. Level spacings s 



m A°AA 



are expressed in units of tA A , and P(s) in units of 

rriAO- AA > W e AA 

FIG. 4. Spectral rigidity: Lines are fits to the form < N 2 > - < N > 2 ~< N > 7 , and the 
values of 7 are 1.036,0.85,0.765 respectively for p* = 0.65, 1.00, 1.35. 

FIG. 5. Participation ratio as a function of unfolded eigenvalues ( (see caption of Fig. 2). 

FIG. 6. The level spacing distribution shown for p* = 1.35 for the full eigenvalue spectrum 
(open rhombs), for the range of eigenvalues with high participation ratio ( 0.75) (filled rhombs), 
along with the RMM prediction. Level spacings s are expressed in units of eA \ — , and P(s) in 

mA(7 AA 

units of 

£AA 
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